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ABSTRACT 


In  the  study  of  some  biological  systems,  control 
theory  may  be  applied.  To  measure  the  transfer  functions  of 
these  biological  transducers,  such  as  a  muscle  spindle,  a  device 
which  produces  a  force  linearly  representing  an  input  electrical 
signal  independent  of  the  compliance  and  damping  of  the  spindle 
is  required.  This  thesis  deals  with  the  design  of  such  a  device 
utilizing  a  Sanborn  pen  motor  as  the  force  producing  device  and 
a  force  transducer.  This  system  was  found  to  be  unstable  due  to 
two  pairs  of  complex  conjugate  poles.  The  system  was  compen¬ 
sated  in  such  a  way  that  the  effects  of  the  complex  poles  were 
almost  cancelled.  Operation  of  the  system  closed  loop  showed  a 
marked  improvement  in  distortion  over  open  loop  operation.  The 
system  was  then  used  to  apply  a  sinusoidal  force  to  a  preparation 
of  frog  muscle  in  an  attempt  to  determine  the  mechanical  small 


signal  transfer  function  of  whole  muscle. 
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CHAPTER  I 


INTRODUCTION 

1.1  PURPOSE 

The  primary  purpose  of  this  project  was  the  building  and 
testing  of  a  control  system  which  would  produce  a  force  output  linearly 
representing  an  electrical  input  signal.  As  the  force  producing  device 
a  Sanborn  type  51-500  galvanometer  was  to  be  utilized.  The  force 
sensing  device  which  would  monitor  the  applied  force  and  also  provide 
a  correcting  signal  was  to  be  a  Statham  type  G  1-4-250  unbonded  strain 
gauge  excited  by  a  Tektronix  type  3C66  carrier  amplifier.  One  app¬ 
lication  of  this  type  of  device  is  in  the  study  of  the  input  versus  out¬ 
put  characteristic  of  mechanoreceptors  such  as  muscle  spindles  and 
cockroach  spines.  The  input  into  such  a  preparation  would  be  the  force 
supplied  by  the  galvanometer;  while  the  output  would  be  the  frequency 
of  nerve  action  potentials.  Since  muscle  spindles  are  said  to  be  non¬ 
linear  with  respect  to  damping  and  spring  constants,  the  secondary 
purpose  of  this  project  was  to  attempt  to  measure  the  passive  trans¬ 
fer  function  of  muscle.  Muscle  was  to  be  tested  rather  than  spindles 
due  to  the  difficulty  of  dissecting  the  spindles. 

1.2  GENERAL 

The  usefulness  of  applying  various  forces  to  biological 
preparations  such  as  muscle  spindles  lies  in  the  fact  that  physiologists 
may  derive  transfer  functions  relating  to  the  individual  preparation; 
such  as  in  Chapman  (3)  and  Chapman  and  Smith  (4).  The  transfer 
function  of  the  me chanoreceptor  studied  helps  in  the  understanding 
of  the  function  of  the  receptor  in  the  maintenance  of  normal  body 


' 


!.  •  r..  - 

. — •  ' 


J. 


jlZQZSl'l, _ I.JL 


M 


t  ",  •,  chi  d  9  *  £".w  *  M.’ '  Q  r  ■  io  :•  ■  •'ll  iq  -.r:  ••■  l 

' 

.  i  '  ' 

, x 

■ 

- 


•;  V  >  . ;  :  V  ■  .  ;  .: >.0  i  ■'  ■■  '  '  *  v  o  £ '4 


,  r- 

. 

L  . 

sJwjy *  sdrxo'tji/I  >T  £  yd  beJbxa  8|j 

-  ff«c  .  /  -  ..'  ■  o  t;b.i  *  «-•  -i-  .  ;.i  *  >/v.r  L»  <»  »  ryj  ax  nox3i»9i<i 


:  -  -[rvao  dbexjra  e;.  vJt,  ao-^OcM.xi  ■>  -•  lo  :  7.,  .  .- XV  i 

>  . 

ov/t  I  '  v  noriJBi  f-  5'-. 'xq  Jn  douz  olni  tvqn:  -at""  ■*e©a;r;a  rfo  r.o 09 

.  i  .  w>V-b  d ;  uri  .  :  9  .-i  v  ■  .  sd  9  -1*  •■ 

' 

Y'x^I>«6o3?  ®xfJ  f  alxtcJando  :>1  iDsaas-*  jrflivr  1& 


,ifl  ;  r  •  ■  »W  .r  s  n  iJ  rat^  '>  ot  Id  .-ton  '  :  v  •'  - 

■  v 

! 

a--  a  n&di  laiil&i  bdtasl  sef  oi  aaw  Maaxx'M  ,  aloaum  lo  xtoi^pcui)  ,f 9I 

1  u"’ 

.  e.*'bn,jqe  &«.xx  >■  ■  ;•  xb  16  ■  xJ  3i!-  •’ 

-  i;  .  .  ;<  id  vol  huos’XMv  gjruyigq.'  -■>  3x;.! :d  s  av  t*dT 

a  ta  “.tcx/vy;  •,.  tsda  1  .  si  «  1  i;i  y*>xl  9.i  -xri  r  ■  y  '  s. f  -?:.•••.  a  ««'  i*-'-  v  0 

!•.  \Uvi.m  ■  re  ixsx/bivib^.r  »di  oJ  5  ■#!  ;o.r izr.  ■  ',0  /.  <  l  .•  r.i 

-  t  V'..  f?  -i  vrfT  ,(^)  ri;r  n3  i£  r  n<?j  3  '  (  .  •-  )-q-.:r) 

:  :5rt&lc  tabiw  9dJ  ni  v.qlsd  ttoJqv  ••  n  »dt  xo  xioiJO'XXjl 

,  r> 

t-. on  :o  »'onto..i9lrtxBrrj  srft  ixi  tol  •-  ’  9  i-  Ic  ioi  ;.  -."1  •' 


2 


function.  The  usual  mechanoreceptors  which  may  be  studied  in  this 
fashion,  such  as  the  spindle  from  the  fourth  toe  of  a  frog's  leg,  are 
the  ones  which  aid  in  reflex  action.  A  classical  example  of  reflex 
action  is  the  simple  knee  jerk. 

The  device  was  designed  to  facilitate  the  application  of 
various  mathematical  forces  to  the  muscle  spindle  which  is  the  sensing 
device  for  the  stretch  reflex  of  the  fourth  toe  of  a  frog's  leg.  Muscle 
spindles  are  bundles  of  modified  muscle  fibres  called  intrafusal  fibres. 
The  sensory  endings  are  coiled  around  a  short  length  of  the  bundle 
and  are  called  anulo  spiral  endings.  The  muscle  spindles  lie  among 
the  ordinary,  or  extrafusal,  muscle  fibres  and  share  their  attach¬ 
ments  so  that  they  are  extended  when  the  muscle  lengthens  and  relaxed 
when  the  muscle  is  shortened.  When  the  muscle  is  elongated,  the 
sensory  portion  is  deformed  resulting  in  depolarization  of  the  excitable 
membranes  and  initiating  a  set  of  action  potentials  .  These  action 
potentials  are  propagated  up  the  group  I  afferent  nerves  coming  from 
the  anulo  spiral  endings.  The  action  potentials  travel  to  a  synapse 
in  the  spinal  cord  and  cause  other  action  potentials  to  be  initiated  and 
propagated  down  the  alpha  efferent  fibres.  These  action  potentials 
reach  the  end-plate  region  of  the  muscle  and  cause  depolarization 
of  the  membrane  which,  in  turn,  causes  the  contraction  of  the 
extrafusal  fibres.  This  type  of  mechanism  has  the  properties  of 
a  length  servomechanism;  for  example,  a  self-regulating  closed- 
loop  mechanism  using  negative  feedback  from  the  spindle  to 
maintain  a  constant  length. 
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A  schematic  diagram  of  a  muscle  spindle  and  its 
attachments  is  shown  in  Figure  1. 


Figure  1.  Muscle  Spindle 

The  mechanoreceptor  part  of  the  above  system  consists  of 
only  the  intrafusal  muscle  fibres,  the  deformation  sensitive  anulo 
spiral  endings  and  the  associated  nerve  fibres.  It  maybe  considered 
as  a  type  of  transducer.  The  transfer  function  of  such  a  transducer 
may  be  derived  in  the  usual  manner  of  applying  a  known  input  and 
studying  the  output  of  the  nerve  fibre;  which  in  this  case  is  impulse 
frequency. 

The  transfer  function  would  be  nerve  impulses  per  second 
divided  by  the  input  force  or  displacement. 

1.3  SYSTEM  REQUIREMENTS 

To  facilitate  the  analysis  of  such  biological  preparations 
the  applied  force  should  linearly  represent  some  input  signal.  This 
means  that  applied  force  should  contain  as  low  a  percentage  of  dis - 
tortion  products  as  possible.  This  requirement  suggests  that  high 

system  gain  and  negative  feedback  must  be  used  to  correct  for  load 
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The  bandwidth  of  such  a  system  should  extend  to  nearly 
one  hundred  hertz  since  prior  work  in  the  field  by  Chapman  and  Smith 
(4)  indicates  that  useful  results  may  be  obtained  in  that  range. 

Useful  force  outputs  would  vary  from  19.  6  x  10”^  newtons 
(20  mg.)  to  about  9.  8  x  10'  newtons. 

The  steady  state  error  for  such  a  system  should  be  as 
s  mall  as  possible  to  eliminate  the  need  for  calibration. 

1.4  DESIGN 

The  usual  method  of  designing  a  servo  system  is  to 
specify  the  required  closed  loop  response  characteristics  and  desired 
accuracy;  then  from  these  requirements  to  choose  components  which 
will  enable  the  system  to  meet  these  specifications.  In  many,  if  not 
all,  cases  some  compromise  must  be  reached  between  specifications 
and  component  costs. 

This  procedure  could  not  be  used  in  this  project  since 
it  was  specified  that  some  of  the  available  components  were  to  be 
used.  These  components  were  the  Gl-4-250  Statham  unbonded  strain 
gauge,  the  Tektronix  carrier  amplifier,  and  the  Sanborn  type  51-500 
galvanometer. 

1.  5  THE  GALVANOMETER 

The  galvanometer  used  in  this  project  was  a  Sanborn 

type  51-500  pen  motor.  It  consists  of  a  permanent  magnet  with  a 

magnetic  field  strength  of  .8  weber/m6.  The  coil  consists  of  3,000 

turns,  is  center  tapped,  and  has  a  resistance  of  3,000  ohms.  It  is 

suspended  in  the  magnetic  field  via  a  torsion  rod  clamped  at  one  end. 

_  3 

The  torque  produced  by  the  galvanometer  is  2  x  10  newton  meters 
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1.6  FORCE  TRANSDUCER 

The  device  which  was  to  be  used  as  the  force  trans¬ 
ducer  is  the  Statham  type  G 1-4-250  unbonded  constantan  strain  gauge 
bridge.  The  force  range  of  this  device  is  1.1  newtons  in  either  direction 
with  a  maximum  excitation  of  9  volts.  Sensitivity  of  the  strain  gauge  is 
25  microvolts /volts  gm.  force,  while  the  compliance  is  .3  mm./kgm. 
force . 

1.7  CARRIER  AMPLIFIER 

To  provide  the  excitation  voltage  to  the  strain  gauge  and 
to  amplify  the  output,  the  Tektronix  type  3C66  carrier  amplifier  was 
used.  This  unit  is  specifically  designed  for  transducer  applications. 

The  excitation  voltage  is  in  the  form  of  a  25  kilohertz  carrier  and  the 
bandwidth  of  this  unit  is  6  kilohertz.  The  amplification  system  is 
essentially  drift-free;  therefore  the  overall  system  drift  is  a  function 
of  the  stability  of  the  external  transducer. 

In  the  operation  of  the  unit,  the  bridge  modulates  the 
applied  25  ikilohertz  carrier  in  accordance  with  the  bridge  unbalance. 
This  bridge  unbalance  is  produced  by  an  input  to  the  transducer  which 
changes  the  resistances  in  the  bridge.  Under  no  input  conditions  to  the 
transducer  the  bridge  is  balanced  and  the  carrier  is  suppressed  at  the 
output.  However,  a  force  input  to  the  transducer  causes  impedance 
changes  which  ,  in  turn,  modulate  the  carrier.  This  a.c.  bridge,  which 
consists  of  the  resistances  of  the  strain  gauge  and  balancing  circuit, 
produces  suppressed  carriers  amplitude  mbdulation. 

The  modulation  side  bands  from  the  bridge  are  applied 
to  an  a.c.  coupled  amplifier  where  the  side  bands  are  amplified 
while  other  frequencies  are  rejected.  These  side  bands  are  then 
applied  to  a  phase  sensitive  demodulator. 
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The  output  of  the  demodulator  is  filtered  and  cor¬ 
responds  to  the  force  signal  applied  to  the  transducer. 

A  block  diagram  of  the  3C66  carrier  amplifier  is 
shown  in  Figure  2. 


(Force)  (Voltage) 


Figure  2.  Block  Diagram  of  3C66  Carrier  Amplifier 
1.8  SPECIFICATIONS 

In  the  design  of  this  system  it  was  expected  that  several 
of  the  requirements  would  be  subject  to  limits  of  equipment.  However, 
several  main  specifications  were  suggested  as  the  extreme  require¬ 
ments  of  the  final  system.  These  specifications  are  as  follows: 

(1)  driving  force  range  9.  8  x  10"^  newtons  (10  mg.  )  to 
9.  8  x  10"^  newtons  (10  gm. )  peak  to  peak  force; 

(2)  bandwidth  D.  C.  to  100  hertz;  and 

(3)  steady  state  error  less  than  2  percent. 
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CHAPTER  II 


2 , 1  OPEN  LOOP  MEASUREMENTS 

In  order  to  obtain  the  characteristics  of  the  system,  it 
was  assembled  as  for  normal  operation  and  the  frequency  response 
was  measured.  To  drive  the  galvanometer  a  driver  amplifier  with 
a  low  output  impedance  was  designed.  The  circuit  was  used  to 
obtain  the  frequency  response  of  the  system  and  it  is  shown  in 
Appendix  I  along  with  gain  and  output  impedance  calculations . 

The  first  frequency  response  curve  was  obtained  with 
a  piece  of  number  18  wire  coupling  the  galvanometer  to  the  transduce 
Other  curves  were  obtained  using  either  the  semi-tendinosus  or 
sartorius  muscle  of  frogs  as  the  coupling  material.  The  galvano¬ 
meter  was  coupled  to  the  transducer  via  the  muscle  in  order  to 
monitor  the  force  applied  to  the  muscle.  In  order  to  prevent 
buckling  the  coupling  material  was  prestretched  by  a  force  of 
about  one  gram  and  the  applied  sinusoidal  force  was  about  9.  8  x  10”^ 
newtons  peak  to  peak. 

The  experimental  arrangements  used  to  obtain  the 
frequency  response  curves  is  shown  in  Figure  3. 


Figure  3.  Block  Diagram  of  Experimental  Set-up  Used  to 
Obtain  the  Frequency  Response  Curves  of  the  System 
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Using  the  experimental  arrangement  shown  in  Figure  3, 
the  open  loop  frequency  response  curves  shown  in  Figure  4  were 
obtained.  The  curves  are  plotted  with  the  gain  normalized  to  one. 

The  shape  of  the  frequency  response  curves  indicates 
that  the  system  consists  of  two  second  order  systems  in  series. 
Additional  information  was  obtained  by  applying  a  step  of  force  to  the 
transducer  showing  a  poorly  damped  system  with  a  resonant  peak  of 
34  decibels  .  The  damping  ratio  calculated  from  the  logarithmic 
decrement  was  found  to  be  approximately  .016.  The  resonant 
frequency  of  this  system  was  measured  to  be  3,000  radians  per 
second. 


The  other  second  order  system  which  occurs  at  the 
lower  frequency  is  much  better  damped.  This  system  consists  of 
the  inertia  of  the  galvanometer,  the  mechanical  characteristics  of 
the  coupling  material  and  the  spring  constant  of  the  transducer. 

Since  the  coupling  material  determines  the  transfer 
function  of  this  part  of  the  system,  the  resonant  frequency  varies 
with  the  coupling  as  shown  in  Figure  4. 

The  differential  equation  which  describes  the  typical 
second  order  system  consisting  of  a  mass  M,  damping  B,  and 
a  spring  constant  K  is  given  as  follows  : 


M  d2x  +  B  dx  +  K  x 
dt  2  dt 


f(t). 


where  f  (t)  is  the  forcing  function. 

Assuming  that  all  of  the  initial  conditions  are  zero 
and  performing  a  Laplace  transformation,  the  following  equations 


result: 


,€  »  x.  j  .  rxi  rwodc  Sn»m*  a£*4&  Is  .  igfcs  t  it  g.'  ' 

■ 

r  %9W  k  •  ni:  nv/odfi  ?9VW3  a*  ••'k,  ot  voitonp  'fi  f'poi  io  OoJ 

,  «»XIO  OJ  r--S;£.[£fr  -lOxI  Cl:  Kg  f)fjj  ,r,f"t'V  ?>:-•*]  oj  ;  €  -  5  c-3  -JD  «••'.  T  A  -  C  •  Cdo 

e^tf,  8 ® v Tii  )■  9  m  xjea*  fp«»trpsii  :  ,rf;s  >, .  ercfsrie  ©ri  [’ 

■  ,  ,  '  <  .  ,  •  r  - 

,sdii»ft  iti  «cnt*>i«Ye  ;©bto  b.  opp  ?  o  J  io  ft.’Riwi*  )  cciaJava  sr  1  isriJ 

S'  i '  ’ . .  ■  ,'j  ’  :  4  -<]'*  yd  bp-M.  -  i'.r-  o  no.-  >  .vrvuVKv 

. 


tc em  :,00  £  >d  0/  b»*ir  v  ni^fc.v?  exifi  lc 

. 


«i5  :  8 1JJD  DO  •.{  >&■£[  ■  '  It;  $3  B\"  •  \*b'lO  fXIO  3  3  IPfJO  »riT 

r  ■'  '  4* 

to  srfaienoo  rnafaye  ?  rH  .•  tcjxusb  *a.tJad  do  -  ?••  yo  a*  r-6>rl  *r o  A'ol  • 

•  i  ■* 

. 

10  v.3ii  Ji?  ji-Tj&dP  isoln.' doe»r>?  oriJ  t‘,53smo£ii3Y.!,  .&g  swr  ;d  :-.  oenron) 

I 

.  : 

•  ■  £  . 

I  ,i  *’1J  Q  3  eSWii/Vtaja  >  ■■'.  r. -••■;■  SI  l  ■ '■' r.  V  ■  P  :  1  ■ 

3.  V")i  -  ri  -  <■£«■  ■•  ’  ,  .pi  ’ c  3c  ■"■  i  '  o-juT.fUo- 

* 

vr  sanal'I  ri.  *u.He  e.&  -  ; ■  'xio  .  -.ri:.-  t  is  •.-■ 

■- 


i>*r.qv;J  £>c?J  a  aob  ri  iri  u-.Ai.  »  ;  i  •  A  •  •»*  oilib  ».f'" 


• 


bn£  ,a  »rttqrrusf>  ,M  t as  x»  e  io  gtriaaleqtto  :n£te  «;i 

1  *  ‘  •  •  .  '  .  '  •  .v-' 

.3  voilo?  qovig  Pi  A  t~uA  no o  £in?q<  t 

,  v  ■  •  ’ 


.0)1  = 


■  ^  r  r,  i 
sb 


ib 


»:o.i Jon 1  r“:  ■  :o-  >o'.  >  .•  -j..  {J?  i  3  vp~: 


4. 


.  -X-  p  j  3(0  •  ii  3  ■  ;i  '•/.  xJ'n  jp  J  j'lb'  I, 


?.->€>;■  up&  gw  wok  o\  oriJ  ,no;.  n  tc ‘.pm  I  3.  1  s  o  to*  •  9'.  p-.n 


*  ■ 


■- 

' 


•.tinea* 


9 


Phase  Shift  (degrees^ 


o 

o 

ro 

i 


Log  Magnitude  db 


10  100  1,000 

Frequency  Hz 

igure  4.  Open  Loop  Frequency  Response  Curves 
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M  s^X  (s)  +  B  s  X  (s )  +  K  X  (s )  =  F(s), 


X  (s) 
F  (s) 


M  s 


+  B  s  +  K  and 


X  (s) 
F  (s) 


K  equals  CO n  which  is  the  natural  undamped  frequency.  B/K  is 

M 

equal  to  2^(Jn.  In  this  case  is  the  ratio  of  the  actual  damping  to 
the  critical  damping  of  the  system.  If  the  equation  is  normalized 
the  following  general  expression  results: 


G  (s)  =  X  (s)  = 
F  (s) 


l*)n 


2  2 

s  +  2  n  +  C*>n 


In  the  case  of  the  transducer, 
COn  =  3,000  radians  per  second. 

G  (s)  _  Eo  /  volti 


=  .016  and 


F: 


in 


Vi 


newt  on 


9.0  x 


9.0  x  10° 

10b  +  97s  +  s2. 


The  other  system  has  a  damping  ratio ^  of.  25  and 
(On  is  1,380  radians  per  second  for  the  wire  coupled  case.  The 
transfer  function  for  this  case  is: 


G(s) 


1.8  x  10 


6 


2  /  6 
s  +  690  s  +1,8  x  10 


'newtons 


volts 


The  total  transfer  function  is  given  by  the  following  equation: 

G  (s)  _  _ 1.8  x  106  x  9.0  x  106 _ 

(s2  +  690  s  +  1.8  x  I06  )  (s2  +  97  s  +  9.0  x  106) 


The  lowest  natural  frequency  measured  was  750 
radians  per  second  using  frog  sartorius  muscle  as  the  coupling 
material.  Changing  the  coupling  resulted  in  very  little  change  in 
damping  ratio:. 
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The  denominators  were  factored  resulting  in  two  pairs 
of  complex  conjugate  poles  for  the  direct  coupled  case  given  below 
as : 

G(s)  _  _ 16.2  x  1012 _ 

“  (s-49.0  +  j3000)  (s-49. 0  -  j3000)  (s-345  +jl295)  (s-345  -  jl295) 

and  for  the  sartorius  coupled  case  as: 

G(s )  _  _ 5.25  x  1012 _  . 

~  (s-49  +  j3000  )  (s -49  -  j3000)  (s  -  187  +j720)(s-187  -j7  20) 

Using  the  above  information  the  root  locus  diagram  shown  in  Figure 

5  was  drawn. 
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igure  5.  Root  Locus  Diagram 
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CHAPTER  III 

DESIGN  OF  COMPENSATION 

3.  1  GENERAL 

The  frequency  response  curve  shown  in  Figure  4 
indicates  that  the  system  is  very  unstable  and  will  not  operate  as 
a  closed  loop  system  except  at  a  very  low  system  gain.  The  static 
loop  sensitivity  calculated  from  the  root  locus  for  stable  operation 
of  the  system  was  found  to  be  K<.  003.  Operation  of  the  system 
closed  loop  at  this  low  gain  offers  no  advantage,  and  in  order  to 
meet  the  specification,  some  method  of  compensation  must  be 
applied. 

In  any  single  second  order  system  containing  complex 
conjugate  poles;  a  lag  network  will  compensate  the  system,  im¬ 
proving  the  steady  state  error  at  the  expense  of  bandwidth.  The 
effect  of  applying  a  lag  to  a  single  second  order  system  is  shown 


System  with  a  Lag  Network 
Since  the  addition  of  a  lag  network  to  an  underdamped 
second  order  system  is  not  adequate,  another  method  of  compensation 
must  be  applied.  This  method  must  in  some  way  cancel  out  the 
effects  of  the  complex  conjugate  poles.  Cancellation  of  these 
unwanted  poles  may  be  carried  out  by  the  addition  of  a  bridged- T 
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network.  This  method  of  compensation  is  described  in  Kuo  (  11) 
and  Chandaket  and  Rosenstein  (  5)  .  The  voltage  transfer  function 
of  a  bridged-T  network  assuming  a  zero  source  impedance  and  an 
infinite  load  impedance  is  given  by  the  following  equation: 

Gc(s)  _  s2  +  2(^/sC4>  +  00 2  . 

s 2  +  2  jsc»>  +  CO  ^ 

A  bridged-T  network  is  shown  in  Figure  7. 
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Figure  7.  Bridged-T  Network 


The  voltage  transfer  function  given  in  terms  of 

circuit  components  is  given  as  follows  : 

E2  (s)  _  1  +  2RC2s  +  R2  C !  C?  s2 
Ei  (s)  “  1  +  R(C!  +  2C2  )s  +  R2  Cl  C2s2  ; 


and 


It  is  seen  from  the  equations  that  the  zeros  may  be 
either  real  or  complex.  The  poles  are  always  real  since  ^>^1 
for  all  values  of  G-  .  For  example,  ^  V  1  =  0  at  C,  =  — L 

^  o  \f  2 

—  yj  2 .  This  places  the  poles  in  a  more  desirable  location  i 


where 


m 


the  s -plane . 


Since  the  system  contains  two  pairs  of  complex 
conjugate  poles,  one  bridged-T  compensator  must  be  designed  so 
as  to  cancel  the  poles  representing  the  force  transducer. 
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Complete  cancellation  is  difficult  and  not  necessary  since  the 
transient  terms  due  to  the  closed  loop  complex  dipole  will  be 
small  if  cancellation  is  incomplete. 

3.2  DESIGN  OF  COMPENSATOR  FOR  THE  FORCE  TRANSDUCER 

To  design  a  compensator  which  cancels  the  unwanted 

poles  of  the  force  transducer,  the  numerator  of  the  compensator  must 

be  equal  to  the  denominator  of  the  transfer  function.  This  sets 

CO  c  =  bJn  and  Cg  c  =  where  C0c  is  the  natural  frequency  and 

c  is  the  damping  ratio  of  the  compensator. 

The  natural  frequency  of  the  force  transducer  is  equal 

to  3,000  radians  per  second  and  the  damping  ratio  is  .016.  Since 

the  compensator  was  designed  for  approximate  cancellation,  (*)c 

was  set  equal  to  3,000  radians  per  second,  and  ^was  set  equal 

to  approximately  .05.  This  results  in  the  following  equations: 

3,000  =  _ 1 _  and  .05  =  /  C^  .  The  solution  of  these 

R  7  C  1C2  Y  Ci 

equations,  which  involves  picking  a  value  for  R  and  solving  for 
Ci  and  C^»  resulted  in  Ci  =  .5  uf,  C 2  =  .001  uf  and  R  =  15  K, 

The  nearest  commercially  available  value  of  Ci  =  .47  uf  was 

i 

chosen.  These  component  values  result  in  c  =  3,  100  radians 

per  second  and  =.046  and  give  the  following  transfer  function: 

Gc  (s)  _  S2  +  285s  +  9.3  x  106 

s2  +  6,700s  +  9. 3  x  106 

The  transfer  function  of  the  bridged-T  network  was 
factored  and  the  root  locus  diagram  was  plotted. 
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The  root  locus  diagram  shown  in  Figure  8  represents 
the  system  compensated  with  a  single  bridged-T  network.  Stability 
of  the  system  was  not  improved.  Static  loop  sensitivity  for  stability 
was  calculated  from  the  root  locus  to  be  K(  3.4  x  10  which  is  the 
value  of  K  on  the  imaginary  axis  of  the  root  locus.  Figure  11  shows 
the  effect  of  the  compensator  on  the  system  in  terms  of  frequency 
response.  For  a  system  gain  of  about  one,  the  system  with  the 
direct  coupling  will  be  stable. 

3.3  DESIGN  OF  A  BRIDGED-T  NETWORK  TO  COMPENSATE  THE 

SECOND  PAIR  OF  COMPLEX  CONJUGATE  POLES 

Since  the  addition  of  a  single  bridged-T  compensator  to 

l 

the  system  showed  insufficient  improvement  in  system  statility, 
another  bridged-T  network  was  added  to  the  system.  This 
compensator  would  have  its  zeros  placed  so  it  would  have  the 
maximum  effect  on  the  two  complex  conjugate  poles  which  are 
dependent  on  the  coupling  material.  A  damping  ratio  for  the 
compensator  of  .5  and  a  natural  frequency  of  1,000  radians  per 
second  were  chosen. 

Solution  of  the  equations,  as  in  the  previous  case, 
resulted  in  the  following  component  values:  C^  —  ♦  1  uf,  Cj  =  .33  uf 
and  R  =5.6  K,.  These  component  values  resulted  in  the  following 
transfer  function  which  is  given  in  factored  form: 

G  r(s)  a  (i  +  540  -  j820)  (s  +  540  +  j820)  . 

(s  +  2,460  )  (  s  +  400) 

The  result  of  this  additional  compensation  is  shown  in 
the  root  locus  diagram  in  Figure  9.  and  in  the  frequency  response 
curve  in  Figure  11.  For  the  system  to  be  stable  the  static  loop 
sensitivity  was  found  to  be  .41.  For  stability  the  system 

gain  from  the  frequency  response  curve  was  about  3.1. 
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igure  8.  Root  Locus  Diagram  Showing  the  System  Compensated 
With  a  Bridged-T  Network 
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igure  9.  Root  Locus  Diagram  Showing  the  Effect  of  Two 
Bridged-T  Compensators  on  the  System 
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3.4  DESIGN  OF  LAG  COMPENSATOR 

For  stability  of  the  system  the  maximum  allowable 
gain  is  about  3.  At  this  gain  the  closed  loop  system  will  have  a 
large  steady  state  error.  To  improve  this,  additional  compensation 
in  the  form  of  a  lag  network  was  applied.  The  two  corner  frequencies 
were  set  at  23  and  330  Hertz  and  the  network  was  incorporated 
in  a  common  emitter  transistor  amplifier  stage.  The  effects  of 
this  compensator  on  the  overall  system  is  shown  in  Figure  10, 
root  locus  diagram,  and  in  Figure  11,  frequency  response  curve. 
Figure  11  shows  the  direct  coupled  case.  Figure  12.  shows  the 
same  compensation  for  the  sartorius  muscle  coupled  case,  as 
well  as  the  closed  loop  response. 

Two  different  muscles  were  used  for  the  closed  and 
open  loop  response  curves.  The  system  was  stable  for  a  gain  of 
about  ten. 

The  transfer  function  of  the  system  is  now  given  by: 

K (s)  Gjq(s)  G2C  (s)  G3<3(s),  where  G  (s)  is  the  original  transfer 
function  and  the  G  nc(s)  are  the  compensators.  For  a  gain  of  12 
the  static  error  coefficient  is  given  as  limit  G(s ) 

A  —  1  U  • 

s - ^0 

This  gain  results  in  a  steady  state  error  of  7 . 7  percent.  A 
complete  evaluation  of  the  system  will  be  made  in  Chapter  6. 

The  design  and  circuit  diagram  of  the  actual  solid  state  circuit 
which  represents  the  amplifier  and  compensators  is  given  in 
Appendix  II. 

3.5  IMPROVEMENT  OF  SYSTEM  RESPONSE 

In  an  attempt  to  improve  the  steady  state  error  the 
corner  frequencies  of  the  lag  network  were  moved  to  lower  fre¬ 
quencies  .  The  transfer  function  of  the  lag  network  is  now 
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igure  10,  Root  Locus  Diagram  Showing  the  Effects  of  Two  Bridged- 

and  a  Lag  Compensator  on  the  System 
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igure  11.  Frequency  Response  Curves  Showing  the  Effect  of 
Compensation  (direct  coupled  system) 
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igure  12,  Frequency  Response  Curves  Showing  the  Open  Loop  and 
Closed  Loop  Response  for  the  Muscle  Coupled  System 


23 


Phase  Shift  (degrees) 
o 


o 

o 

t-H 

I 


o 

o 


o 

o 

o 


o 

o 

s^H 


o 

s=H 


si 

X 

u 

d 

CD 

2 

cr 

D 


Ll 


igure  13.  Frequency  Response  Curves  for  the  Open  Loop  and 
Closed  Loop  System  with  the  Corner  Frequency  of 
the  Lag  Network  Moved  to  ]  0  Hz 
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Gc(s)  =  (s  +  1,000)  . 

(s  +  63 ) 

The  effect  of  this  change  is  shown  in  the  frequency  response  curve 
in  Figure  13.  For  the  direct  coupled  case  at  a  gain  of  14  the  phase 
margin  is  measured  to  be  about  40°. 

3.6  MINIMIZATION  OF  DISTURBANCES 

In  the  operation  of  the  system  at  low  force  levels,  in  the 
order  of  30  to  50  mg.  ,it  was  noted  that  excessive  noise  existed  making 
work  at  these  levels  useless.  The  majority  of  the  noise  came  from 

the  3C66  carrier  amplifier  due  to  imperfect  filtering  of  the  carrier 

* 

and  harmonics.  The  output  level  of  noise  did  not  depend  on  the 
carrier  amplifier  gain.  Another  form  of  noise  came  from  building 
vibration. 

To  minimize  the  effects  of  these  two  sources  of 
extraneous  disturbances,  a  method  of  gain  distribution  was  used. 

A  block  diagram  showing  the  disturbances  with  respect  to  the 
system  is  shown  in  Figure  14. 


Figure  14.  A  Block  Diagram  Showing  the  Disturbances 

Affecting  the  System 

For  the  first  case  consider  that  if  N -p  =  N?  =  0,  then 

F(s )  =  AG  (s)  E  (s)  . 

1  +  AG  (s) 
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If  E(s  )  —  0,  N,  J.  0,  and  N2  =0;  then  F(s)=  NA _ 

1  +  AG(s ) 

Therefore,  in  order  to  minimize  the  effects  of  Nj,  the  ratio, 

Nl  ,  must  be  made  as  small  as  possible. 

G(s )  E(s ) 

To  minimize  the  effects  of  N2  the  ratio  of  _ £i2 _ 

AG(s )  E(s ) 

must  be  minimized.  In  the  actual  system  this  was  achieved  by  having 
a  high  gain  in  the  forward  amplifier  stage.  The  gain  of  the  carrier 
amplifier  was  made  very  high  and  the  output  attenuated  so  that 
the  disturbance,  which  was  independent  of  gain,  was  attenuated 
by  some  factor. 

In  the  final  design  a  selector  switch  was  inserted  at 
the  output  so  that  the  gain  of  the  carrier  amplifier  could  be  main¬ 
tained  at  the  required  levels  for  total  system  gain  by  the  adjustment 
of  carrier  amplifier  gain  and  attenuation. 
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CHAPTER  IV 

MEASUREMENT  OF  MUSCLE  TRANSFER  FUNCTION 

4.1  OBJECT 

The  secondary  objective  of  this  project  was  to  use  the 
servo  system,  as  designed  in  the  first  part,  to  define  the  biological 
load;  i.e.,  to  measure  the  passive  mechanical  transfer  function 
of  skeletal  muscle.  This  study  was  carried  out  on  the  striated 
or  skeletal  muscle  of  frogs. 

4.2  MUSCLE  PHYSIOLOGY 

Striated  muscle  consists  of  individual  fibres  which  are 
ten  to  one  hundred  microns  in  diameter.  The  individual  fibres 
are  surrounded  by  a  polarized  membrane,  the  inside  of  which  is 
about  one  hundred  millivolts  negative  with  respect  to  the  outside. 

Each  fibre  is  made  up  of  a  number  of  parallel  elements 
about  one  micron  in  diameter  which  are  called  myofibrils  .  The 
myofibrils  themselves,  are  made  up  of  parallel  actin  and  myosin 
filaments  in  a  regular  array.  The  myosin  filaments  are  about  one 
hundred  and  sixty  angstroms  in  diameter  and  about  one  or  two 
microns  in  length.  The  actin  filaments  are  about  one  micron  in 

I 

length  and  fifty  to  seventy  angstroms  in  diameter.  The  overlap 
between  the  arrays  of  actin  and  myosin  give  rise  to  the  regular 
pattern  of  striations  as  viewed  under  a  microscope.  A  schematic 
illustration  of  a  myofibril  is  shown  in  Figure  15. 
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Actin  and  Myosin  Filaments  in  a  Myofibril.  (  The 
thick  filaments  are  myosin  and  the  thin  filaments 
are  actin). 

The  pattern  shown  in  Figure  15  is  characterized  by  a 
succession  of  bands  (dense  bands  called  A  bands  and  light  bands 
called  I  bands).  The  z  line  is  a  narrow  structure  to  which  the  actin 
filaments  are  attached.  It  has  been  reported  by  H.E.  Huxley  (7) 

(  8  )  that  cross  bridges  exist  between  the  thick  and  thin  filaments 
which  are  the  only  mechanical  linkage  between  the  filaments. 

Contraction  of  muscle  is  preceded  by  depolarization 
of  the  membrane.  This  depolarization  causes  the  release  throughout 
the  fibre  of  an  activating  substance  which  is  now  accepted  to  be 
calcium  ions.  The  release  of  calcium  enables  one  of  the  proteins, 
myosin,  to  act  as  an  enzyme  and  split  a  phosphate  group  from 
adenosine  triphosphate  (ATP)  which  provides  the  energy  for 
contraction.  Adenosine  triphosphate  is  a  high  energy  active 
phosphate  compound  which,  by  losing  one  of  its  phosphate  groups, 
becomes  adenosine  diphosphate  (ADP)  to  supply  energy  to  drive 
chemical  reactions  which  otherwise  would  not  proceed. 
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The  complete  mechanism  of  contraction  is  by  no  means 
clear,  but  it  is  now  thought  that  the  actin  and  myosin  filaments  slide 
past  each  other  producing  changes  in  length  of  the  fibre.  It  has 
been  suggested, and  to  some  extent  verified  by  H.E.  Huxley  (7)  (8) 
that  a  relative  force  is  exerted  at  the  sites  of  the  cross  bridges 
resulting  in  this  sliding  action. 

Mechanically,  the  result  of  a  striated  muscle  being 
stimulated  is  a  "twitch",  which  is  characterized  by  a  phase  of 
contraction,  followed  by  a  phase  of  relaxation.  If  the  muscle  is 
held  at  a  fixed  length  a  tension  change  will  be  recorded  during  the 
twitch.  If  the  muscle  is  stimulated  again  before  it  has  completely 
relaxed  the  second  contraction  will  be  superimposed  on  the  first 
contraction.  Rapid  stimulation  which  does  not  allow  the  muscle 
to  relax  results  in  prolonged  mechanical  response  called  a 
tetanus . 

It  is  generally  accepted  (  10  )  (  12)  that  active  muscle 
consists  of  a  passive  elastic  component  and  a  contractile  component 
in  series.  The  contractile  component  consists  of  a  parallel 
combination  of  damping  and  elastic  components.  The  series 
elastic  and  parallel  damping  component  tend  to  distort  the 
rapid  changes  in  tension. 

Two  possible  mechanical  models  of  muscle  in  the 
active  state  have  been  proposed  by  Jewell  and  Wilkie  (10). 

In  the  unexcited  state  muscle  also  contains  the  series 
elastic  component  (12).  The  parallel  elasticity  and  damping  exist 
but  probably  bear  little  or  no  relation  to  these  components  in  the 
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active  state  (2).  In  whole  muscle  the  membranes  would  tend  to 
contribute  to  the  damping  and  compliance. 

4.  3  GENERAL 

Since  a  muscle  is  part  of  the  overall  control  system, 
it  would  be  useful  to  describe  it  in  terms  of  a  mechanical  transfer 
function.  The  knowledge  of  the  transfer  function  of  the  coupling 
material  would  aid  in  the  determination  of  overall  system  stab¬ 
ility  and  it  would  contribute  some  information  about  the  proper¬ 
ties  of  passive  muscle. 

Although  skeletal  muscle  is  studied  in  this  case, 
it  is  not  necessarily  mechanically  equivalent  to  the  muscle 
spindle.  Skeletal  muscle,  however,  is  the  best  alternative  to 
the  muscle  spindle  and  much  easier  to  work  with  and  dissect  out. 

4.4  PROCEDURE 

The  usual  method  of  analysis  of  a  mechanical  system 

consisting  of  mass,  damping  and  compliance  is  to  state  the 

displacement  at  some  point  as  the  output  with  force  as  the  input. 

The  transfer  function  is  given  as  where  X(s)  and  F(s)  are 

F(s) 

the  Laplace  transforms  of  the  displacement  and  force  respectively. 
This  relation,  of  course,  assumes  a  linear  transfer  function. 
Previous  work  by  Buchthal  (2)  indicates  that  a  linear  relation  may 
be  obtained  for  single  and  small  bundles  of  fibres  if  small  variations 
of  force  are  used.  From  this  it  is  clear  that  in  order  to  obtain  a 
transfer  function  of  whole  muscle,  a  known  force  must  be  applied 
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to  a  muscle  and  the  displacement  measured. 

4.  5  FORCE  MEASUREMENT 

The  Statham  Gl -4-250  force  transducer  was  used  to 
monitor  the  force  applied  to  the  muscle.  If  the  mass  of  the  muscle 
is  small,  an  equal  and  opposite  force  appears  at  the  transducer, 
monitoring  the  force  applied  by  the  galvanometer.  In  this  way,  the 
force  applied  by  the  galvanometer  may  be  monitored  until  the 
force  due  to  the  mass  and  acceleration  of  the  muscle,  introduces 
a  phase  lag  between  the  applied  and  measured  force.  This  effect 
depends  on  the  effective  mass  and  the  compliance  of  muscle;  thus 
if  the  measurements  are  made  at  frequencies  below  the  resonant 
frequency  of  the  muscle,  little  or  no  error  is  introduced  in  the 
measurement  of  the  force, 

4.  6  DISPLACEMENT  MEASUREMENT 

The  displacements,  or  length  changes,  of  the 
muscle  were  measured  using  a  light  beam  lever  as  shown  in 
Figure  16, 


Figure  16.  Displacement  Measurement  System 
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The  light  source  was  focused  on  a  light  mirror  attach¬ 
ed  to  the  galvanometer  pen  and  reflected  on  to  a  solar  cell,  so  that 
changes  in  the  length  of  the  muscle  resulted  in  changes  in  the  area 
of  the  solar  cell  illuminated.  The  rise  time  of  the  solar  cell  is 
given  as  20  microseconds.  The  system  was  calibrated  against  a 
dissecting  microscope  which  measured  the  displacement  directly. 
4.7  PROCEDURE 

The  muscle  was  dissected  from  freshly  killed  frogs 
and  the  tendons  were  tied  with  steel  wire  and  attached  via  this 
wire  to  the  galvanometer  and  the  force  transducer.  The  sartor  - 
ius  and  semi-tendinosus  muscles  of  frogs  were  used.  In  one 
case  the  semi-tendinosus  muscle  of  a  toad  was  used.  During  the 
test  the  preparation  was  kept  moist  by  continual  application  of 
Ringer’s  solution.  The  measurements  were  made  on  muscles 
weighing  250  mg.  or  less.  The  toad  semi-tendinosus  muscle 
weighed  250  mg.  and  had  an  estimated  resonant  frequency 
of  about  200  hertz . 

Using  the  experimental  arrangement  described, 

the  frequency  response  was  measured  over  a  frequency  range 

*• 

from  .  01  to  100  hertz.  The  system  was  operated  closed  loop  and 
'  exhibited  the  closed  loop  frequency  response  shown  in  Figure  13 
(wire  coupled  case).  The  measurements  were  independent  of 
system  characteristics  due  to  the  fact  that  both  the  output  and 
input  were  measured.  The  magnitude  of  the  force  applied  was  kept 
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constant  over  the  range  of  measurement  to  eliminate  errors  due 
to  the  nonlinearity  of  compliance  of  the  muscle.  This  was  done 
by  varying  the  input  voltage  which  provided  the  reference  input 
to  the  control  system.  The  relative  magnitudes  of  input  (refer¬ 
ence),  force  out,  displacement  and  the  phase  shifts  of  the  outputs 
were  measured  as  in  Chapter  II. 


A  diagram  of  the  experimental  arrangement  is 
shown  in  Figure  17. 


Figure  17.  Block  Diagram  of  Experimental  Arrangement 

Used  to  Measure  the  Mechanical  Characteristics 

of  Muscle 

Using  the  arrangement  shown  in  Figure  17, 
frequency  response  curves  were  obtained  for  several  muscles. 

Some  initial  tension  was  applied  to  the  muscle  and  the  magnitude  of 
the  force  variation  was  kept  small  so  that  the  nonlinear  effects 
of  muscle  compliance  would  have  little  effect  on  the  variation 
(i.  e.  ,  the  compliance  would  have  essentially  a  constant  value  over 
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the  complete  variation).  The  nonlinear  characteristic  of  muscle 
compliance  is  discribed  by  Buchthal  (2)  and  Wilkie  (12). 

4.8  RESULTS 

The  measurements  on  muscle  yielded  the  magnitude 
and  phase  shift  of  the  applied  force  and  the  resulting  displacement 
with  respect  to  the  reference  voltage  (input). 

Normalized  curves  of  are  plotted  in  Figures 

18  and  19.  Curve  1  was  obtained  using  frog's  sartorius  muscle 
which  was  prestretched  to  a  test  length  of  1.13  times  the  length 
in  the  body;  while  in  curve  2  the  same  muscle  was  prestretched 
to  1.  35  times  the  initial  length.  The  weight  of  this  sample  was 
150  mg.  Curves  3  and  4  were  obtained  using  frog's  semi- 
tendinosus  muscles  prestretched  to  a  length  of  about  1.  1  times 
the  initial  length  and  weighing  100  and  94  mg.  respectively. 

Curve  5  was  obtained  using  a  toad's  semi-tendinosus  muscle 
weighing  250  mg.  which  was  prestretched  to  a  length  of  1 .  1 
times  the  initial  length.  The  estimated  resonant  frequency  of 
this  sample  was  about  200  hertz.  Curve  6  was  obtained  using 
frog's  semi-tendinosus  muscle  weighing  110  mg.  and  prestretch¬ 
ed  to  a  length  of  1.15  times  its  initial  length.  The  above  results 

_  o 

were  obtained  using  a  sinusoidal  force  less  than  9.  8  x  10 
newtons  peak  to  peak. 
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igure  19.  Frequency  Response  Curves  of  Typical  Muscle 
Transfer  Function  Measurement  Results 
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4.  9  DISCUSSION 

The  results  obtained  from  the  measurement  of  the 
mechanical  transfer  function  of  muscle  shown  in  Figures  18  and 
19  do  not  represent  any  simple  linear  transfer  function.  These 
results  do  not  compare  with  the  linear  relationships  obtained  by 
Buchthal  (2)  on  bundles  and  single  fibres.  The  results  obtained 
by  Buchthal  (2)  show  clearly  one  corner  frequency  which  would 
be  due  to  a  parallel  combination  of  compliance  and  damping.  No 
such  distinct  corner  is  evident  in  the  results  shown  in  Figures 
18  and  19.  Howe\er  at  certain  frequencies  a  sudden  change  in 
phase  is  evident  in  nearly  all  the  tests  done.  Buchthal  (2)  how¬ 
ever,  shows  no  measurements  of  phase  shift. 

One  reason  for  these  unusual  results  could  be  that 
compliance  of  the  connective  tissue  masked  the  compliance  and 
damping  of  the  muscle  fibres.  This  does  not  explain  the  sudden 
change  in  phase. 

Another  possibility  which  seems  quite  unlikely  is 
that  the  mechanical  oscillation  causes  excitation  of  the  fibres  and 
hence,  contraction.  An  attempt  to  eliminate  this  possibility  by 
prestretching  the  muscle  so  that  the  actin  and  myosin  filaments 
no  longer  overlap  failed  due  to  the  large  forces  needed  to  accomplish 
this.  Under  these  prestretched  conditions  the  muscle  would  no 
longer  contract.  I  f  excitation  did  indeed  exist ,  tetanus  might  ex¬ 
plain  the  change  in  phase;  however,  the  results  do  not  confirm  this. 
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CHAPTER  V 

EVALUATION  OF  SYSTEM 

5.1  DIRECT  COUPLED  CASE 

A.  Transient  Response 

A  5  hertz  square  wave  was  applied  to  the  compensated 
system  with  the  open  loop  gain  set  at  10.  The  system  was  direct 
coupled  and  compensated  as  shown  in  the  frequency  response 
curves  in  Figure  11.  The  square  wave  response  is  shown  in 
Figure  20.  The  damping  ratio  for  this  gain  is  calculated  to  be  about 
.45. 


Figure  20.  Transient  Response  of  Direct 
Coupled  System 

B.  Steady  State  Error 

With  the  system  gain  equal  to  10  the  static  error 
coefficient  of  the  type  zero  system  K0  is  defined  as: 


Kq  =  limit.  K  G(s )  Gcl(s)  G  2(s)  Gc3  (s). 
s — >-0 

For  the  system  in  question,  K0  —  10  since  K  =  10. 
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The  steady  state  error  is  given  as  follows: 


ess 


_ R _ 

1  +  Limit  G(s)  <Gc(s) 

s — >0 


where  egs  =  100% 

1  +  K0 

=  9%. 


The  steady  state  error  measured  from  the  square  wave  input  was 
also  about  9  percent.  The  gain  was  increased  to  15  and  the 
steady  state  error,  both  measured  and  calculated,  was  about 

6  percent.  This  mode  of  operation  was  poor  for  frequencies  over 
30  hertz  since  the  damping  ratio  was  very  low  resulting  in  a  larg 
resonant  peak. 

C .  Bandwidth 

The  bandwidth  of  the  system  with  a  gain  of  10  is 

7  0  hertz.  A  resonant  peak  is  evident.  The  bandwidth  may  be 
increased  by  an  increase  in  gain  with  an  associated  increase  in 
the  resonant  peak. 

D.  Distortion 

Distortion  measurements  on  the  system  were  made 
at  20  hertz  and  50  hertz.  The  distortion  measurements  are  giyen 
below  in  Table  I. 


F  requency 

Open  Loop 
Distortion  % 

Closed  Loop 
Distortion  % 

20  hertz 

3.5 

.7 

50  hertz 

4.5 

1.8 

Table  I.  Comparison  of  Closed  and  Open 
Loop  Distortion. 


E.  Drift 

In  the  design  of  the  circuit,  silicon  transistors 


were  used  and  the  stages  were  alternate  pnp  and  npn  to  minimize 
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temperature  drift.  The  force  transducer  was  used  in  the  four 
arm  bridge  configuration  in  order  to  minimize  temperature  effects. 
No  measurements  of  temperature  sensitivity  were  made  because 
the  final  system  was  to  operate  closed  loop,  and  the  carrier 
amplifier  and  transducer,  which  contribute  to  the  overall  drift, 
were  purchased  items. 

5.2  MUSCLE  COUPLED  CASE 

A.  Transient  Response 

The  effect  of  a  5  hertz  square  wave  on  the  system 
with  a  gain  of  10  is  shown  in  Figure  21. 


Figure  21.  Transient  Response  of  Muscle  Coupled  Cases 

The  damping  ratio  was  calculated  to  be  .34.. 

B.  Steady  State  Error 

The  steady  state  errof  to  a  step  is  the  same  as  in 
the  direct  coupled  case  which  is  9  percent. 

C.  Bandwidth 

The  -3  db  point  with  a  gain  of  10  was  65  hertz  and 
the  resonant  peak  was  measured  to  be  5  db  at  50  hertz.  The  results 
of  these  measurements  are  given  in  Table  II. 
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D.  Distortion 

Distortion  measurements  were  again  made  on  the 
system  at  20  hertz  and  50  hertz. 


Frequency  Open  Loop  Closed  Loop 

'  ”  Distortion'%  Distortion  % 

20  hertz  6.5  1 

50  hertz  7  4.7 

Table  II.  Comparison  of  Closed  and  Open  Loop  Distortion. 

5.3  LAG  NETWORK  AT  LOWER  FREQUENCY 

A.  Transient  Response 

For  a  gain  of  14,  the  damping  ratio  was  calculated 
to  be  .50.  The  step  response  is  shown  in  Figure  22. 


Figure  22.  Transient  Response  with  Lag  Network 
at  Lower  Frequency 

B.  Steady  State  Error 

The  steady  state  error  for  this  system  was  measured 


to  be  about  6  percent. 
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C.  Bandwidth 

The  bandwidth  for  this  case  was  about  7  0  hertz. 

The  resonant  peak  was  about  4  db. 

D.  Distortion 

The  open  loop  distortion  is  5  percent  as  compared  to  .9 
percent  closed  loop  at  20  hertz. 

5.4  CONCLUSION 

The  original  system  which  was  very  unstable  was 
compensated  using  two  bridged- T  networks  and  one  lag  network. 

The  effect  of  this  compensation  resulted  in  the  steady  state  error 
of  6  to  7  percent.  This  did  not  meet  the  original  specification  of 
2  percent.  If  a  low  steady  state  error  is  required,  the  lag  network 
could  be  set  to  a  lower  frequency  resulting  in  a  higher  gain  but 
with  a  resulting  decrease  in  bandwidth.  The  damping  ratio  likewise 
depends  on  gain  and,  if  the  system  is  to  be  used  with  a  step  function 
as  the  input,  the  gain  could  be  lowered  to  decrease  the  bandwidth 
and  to  increase  the  steady  state  error  with  an  improvement  in  the 
damping  ratio. 

One  of  the  main  reasons  for  the  design  of  this  system 
was  to  decrease  the  distortion  of  the  applied  force  due  to  muscle 
nonlinearity.  The  amount  of  improvement  is  dependent  on  the  gain 
and,  from  the  results  of  measured  distortion,  a  definite  improvement 


is  seen. 


•i-idri  OT  JuocU  3,53  ftirii  10I  dibivfbnsd  odl 

.db  £*  3;f'jdfc.  ?.f.  v  >  c  •'  •;■  i<  '  d^O^T 


1 

toriro  ?.eK7  -C 


r  ,s  •.  u  a  a  2ivk>  ■  ■•  ’  ai  tc  •  ;X  ,J..  C;^>o:  t  •  '  i 


j 


s  ;  ad  OS  u  -  >oJ  baaoja  tnoDtoq 

^lftV,42HQ2 - £*& 

g.  1W  -Icfr^afUr  ,  T-  *  ::*iiiw  rrf»38>fe  Isdigio  >nT 

-to  CTO  O)  >  63JP.  e*:i  m  Mth/W*  a  :i£L  ?i}  eirii  1©  t©  ^  9rfl 

*  ^ 

. 

i 

tr;d  nifSS  T^irf  £ 

" 

- 

Qtii  d,  j«dmavo7qttix  as  dtiw  lotto  •#■  »  9dt  ooso'ton:  c*  b-n* 

.  '  *  _  *Ja  1  -  ■'  V 


alDsara  o3  90b  9  to}  b/uiqqs  arf*  \o  roiJioJatb  «dl  9WMWb  o*  •£’* 

■ 

rtarnavo-rqmi  at ;n  i  >b  '  b9-M**mio  *3b.9v?  r  rto'a  «'  a~ 


.ns'jfe  ai 


42 


For  use,  at  up  to  30  or  40  hertz,  with  a  sinusoidal 
input  the  system  may  be  operated  at  higher  gains.  (A  sinusoidal 
input  may  be  used  in  the  study  of  mechanor eceptor s .  ) 

5.  5  RECOMMENDATIONS 

As  suggested  in  Chapter  I,  it  was  found  that  the 
equipment  limited  the  achievement  of  the  desired  specifications. 
Since  the  system  does  not  quite  meet  all  specifications  outlined 
in  Chapter  I,  some  suggested  improvements  could  be  the  use  of 
a  faster  galvanometer  and  as  the  force  transducer,  the  use  of  a 
unit  with  a  much  higher  resonant  frequency  and  a  better  damping 
ratio. 
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APPENDIX  I 

In  order  to  obtain  the  frequency  response  curves  of  the 
system,  an  amplifier  was  built  to  drive  the  galvanometer.  The  output 
stage  consisted  of  a  Darlington  arrangement  used  as  an  emitter 
follower.  The  input  stage  was  a  differential  amplifier.  A 
circuit  diagram  is  shown  in  Figure  A-l. 
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Ti,  T2  ,  T3  ,  T4 
T5  =  TI  487 
R1  =  R2  =  33K 
R3  =  2K 
Rc  =  1.5K 


TI  4  95 


R7  =  5 . 6K 
R6  =  12K 
R5  =  12K 


Driver  Circuit  Used  to  Obtain  the 
Frequency  Response  Curves  of  the 
System 


Figure  A  - 1 , 
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The  output  impedance  of  the  configuration  is  given  by 


JLLl 


hfe4  ^fe5 


the  expression  Z  out  ^  ^ib5  + 

26mv  +  1.5  ^  2.8  ohms. 

Therefore  Z  out  ^  10  ohms. 


where  h 


ibc  = 


JE5 


The  single  ended  gain  of  the  circuit  is  given  by  the 


following  expression: 


Av  = 


R 


hib  i  +  hib2  +  rE 
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APPENDIX  II 

The  compensators  as  designed  in  Chapter  3,  and  the 
driver  amplifier  shown  in  Appendix  I,  were  incorporated  into  a  single 
amplifier.  Also  included  in  the  circuit  was  a  differential  amplifier 
needed  to  provide  the  error  signal. 

INPUT  STAGE 

To  provide  an  error  signal  to  drive  the  galvanometer  a 
differential  amplifier  was  used.  It  was  designed  with  a  temperature 
compensated  current  source  and  emitter  resistors  to  aid  in  the  D.C. 
balance  of  the  system.  The  input  stage  is  shown  in  Figure  A -2. 


The  voltage  gain  for  a  single  ended  output  is  given  by: 
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In  the  actual  circuit  the  voltage  gain  is  approximately 


1. 2  for  an  operating  current  of  .  5  ma. 


The  single  ended  input  impedance  is  given  by: 


Zin  =  hfe  (2hib  +  Re)  // 


4  megohms. 


IMPEDANCE  MATCHING  OF  COMPENSATION 


In  order  not  to  affect  the  transfer  function  of  the 


bridged- T  compensators  the  input  and  output  impedances  were 
matched  using  emitter  followers.  Alternate  pnp  and  npn  stages 
were  used  to  minimize  temperature  effects. 

LAG  NETWORK 

The  lag  network  was  incorporated  into  a  common 
emitter  transistor  stage  as  shown  in  Figure  A-3. 


Figure  A-3. 
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Gain  of  the  stage  is  2.  2  at  D.  C. 

OPERATIONAL  AMPLIFIER  STAGE 

In  order  to  provide  additional  gain  to  the  system 
which  could  be  adjusted  without  shifting  the  D.  C.  level  a  simple 
operational  amplifier  was  used.  A  circuit  diagram  is  shown  in 
Figure  A-4. 


R 


L2 


Figure  A-4.  Operational  Amplifier 

The  operating  current  of  the  differential  stage  was 

set  at  .  5  ma. 


The  open  loop  gain  of  this  amplifier  is  given  as 


follows: 


Av  = 


R 


x 


L2  ^  450. 


2h 


+  R 


R. 


ib  *'  E  "*E 

The  diode  was  again  used  to  temperature  compensate 

the  current  source. 

This  amplifier  was  then  followed  by  the  driver  stage 
which  is  discussed  in  Appendix  I.  A  complete  circuit  diagram  is 
given  in  Figure  A-5. 
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igure  A-5.  Circuit  Diagram  of  Complete  Circuit 
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LIST  OF  COMPONENTS 

T  ,  T2,  T?,  Tg  =  2N3707 
Tg,  T5,  T9  =  2N3705 
T4>  T6,  T1q  =  2N3703 


Tn  =  TI415 

T12  =  TI487 

Dlf  D2  =  1N1020 

Rx  =  15  M 

R9  =  15  K 

Rj7  =  12  K 

R2  =  27  K 

R10  =  39  K 

R18  =  12K 

Rg  =  2.7  K 

Rn  =  5.  6  K 

R19  =  5.6  K 

R4  =  20  K 

R12  =  6.  8  K 

R2q  =  10  K 

R5  =  12  K 

R 1 3  =  15  K 

R21  =1M 

R6  =  12  K 

R14  =  1  K 

R22  =  2.  2  k 

R?  =  5.  6  K 

R15  =  2  K 

R23  =  33  K 

Rg  =  15  K 

Cj  =  .  001  uf 

C2  =  .  5  uf 

Cg  =  .1  uf 

C  .  =  .  33  uf 

4 

Cg  =  .68  uf 

+Vg  =  20  v 

-Va  =  20  v 
s 

Rl6  =  4.  7  K 

R  =  1.  5  K 
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